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Abstract 


This article describes the fluid dynamics video: "Effect of bubble deformability on 
the vertical channel bubbly flow." The effect of bubble deformability on the flow rate 
of bubbly upflow in a turbulent vertical channel is examined using direct numerical 
simulations. A series of simulations with bubbles of decreasing deformability reveals 
a sharp transition from a flow with deformable bubbles uniformly distributed in the 
middle of the channel to a flow with nearly spherical bubbles with a wall-peak bubble 
■ distribution and a much lower flow rate. 

1 Introduction 

A series of numerical simulations of an upward bubbly flow with a void fraction of 3% 
between two vertical walls using direct numerical simulation with a front tracking method 
is performed. A rectangular computational domain with two no-slip vertical wall and pe- 
riodic boundary conditions in streamwise and spanwise directions is used. The effect of 
bubble deformability is isolated through controlling the surface tension while keeping other 
flow parameters fixed. The flow is driven upward by a specified pressure gradient. The 
dimensionless numbers in the problem are Reynolds number, Re = E^HL^ Eotvos number, 

Eo = and density and viscosity ratios. Here, p, [i, a, d, w, u, and g are the liquid density 
and viscosity, surface tension coefficient, bubble diameter, channel width, average vertical 
velocity in the channel, and acceleration due to gravity, respectively. The pressure gradient 
and Eotvos number are specified a priori and the Reynolds number is calculated from the 
average velocity a posteriori. Eotvos number is decreased from 4.5 to 0.9. The decrease 
in the deformability of bubbles changes the lateral lift force on the bubbles. Therefore, 
the distribution of the bubbles changes from a uniform distribution in the middle of the 
channel to a wall-peak distribution with only a few bubbles in the center. The reduction of 
the number of the bubbles in the center increases the average density there which leads to 
a reduction in flow rate since the driving pressure gradient is kept fixed. As a result, the 
Reynolds number drops from 3750 to 1700. 
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